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COMPUTER SIMULATION INCORPORATING A HELICOPTER MODEL FOR 


EVALUATION OF AIRCRAFT AVIONICS SYSTEMS 

Aaron J. Ostroff and R. Brian Wood* 
Langley Research Center 


SUMMARY 


A computer program has been developed to integrate avionics research in nav- 
igation, guidance, controls, and displays with a realistic aircraft model. A 
user oriented program is described that allows a flexible combination of user 
supplied models to perform research in any avionics area. A preprocessor tech- 
nique for selecting various models without significantly changing the memory 
storage is included. Also included are mathematical models for several avionics 
error models and for the CH-47 helicopter used in this program. 


INTRODUCTION 

Avionics research is being conducted at the Langley Research Center in 
areas related to navigation, guidance, control, and displays (ref. 1). An over- 
all simulation is required to evaluate and integrate the various areas of this 
research. The purpose of the present paper is to describe the simulation pro- 
gram and its capabilities. 

One requirement is to develop a computer program that gives each researcher 
the capability of evaluating an avionics area independently from other research 
areas without rigid time schedules. Further, algorithms that are developed 
should be included into the simulation with minimal impact on the overall pro- 
gram. The preceding requirement has been met by using a modular type design. 

A set configuration representing all of the avionics systems envisioned has been 
defined, with separate modules representing the various measurement sensors and 
avionics functions. Potentially, 10 different user supplied models can be 
selected for each of the modules in the simulation. The choice of models is 
accomplished by a series of defaults; overriding can be accomplished with the 
use of input data cards. A specially designed preprocessor combines the input 
data and default data in order to assemble a complete program. This program 
includes a six-degree-of-freedom CH-47B helicopter model (ref. 2). After a 
model has been developed, it is added to the program subroutine library for use 
by other researchers. 

The computer program must be simple to use and easily understood since 
users of various disciplines will be running the program. NAMELIST type data 
are used for the input, and only the data associated with the particular simula- 
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tion being run are required. This feature is accomplished by the series of 
defaults previously described. 

Another requirement that the computer program should have is the capability 
of replacing any computer software module by its equivalent hardware and of eval- 
uating the hardware in the overall simulation system. In certain cases, pilot 
capability is also required. These features dictate the need for a real-time 
simulation. Two programs have been developed and are maintained, one for batch 
type research and the other for real-time capability. These programs are main- 
tained as nearly similar as possible within the limits of the Langley computer 
complex. The real-time program can be combined with hardware to meet this 
requirement. The computer software module can be deleted by a simple input data 
card that selects a dummy model with zero iteration frequency. The equivalent 
input communication is carried out through electrical input/output channels to 
the actual hardware. 

Elements of structured type programing have been used where applicable. 

The mainline and much of the programing in the subroutines use top-down struc- 
turing. With the Langley FORTRAN compilers, structured programing tends to be 
inefficient with respect to time and, to some degree, memory storage. The effi- 
ciency problem has been traded off to maintain single- input, single-output type 
programing as much as possible. 

An overlay program has been developed to minimize core storage. In gen- 
eral, the program can be run with approximately 60 000 octal words. The itera- 
tion frequency of every module is selectable by the user through a simple data 
card. Print options are also selected by input data cards, with each module hav- 
ing at least one option. 


SYMBOLS 


Ai,A 2 ,A 3 ,Ai| matrix components representing stability derivatives (eq. (15)) 
B-|,B 2 ,B 3 ,B 4 matrix components representing control derivatives (eq. (15)) 

Cg matrix transformation from body frame to geographic frame 

C b matrix transformation from geographic frame to body frame 


d^ correlation distance, m (eq. (12)) 

E error term 


F 9 * 9 system matrix (eq. (10)) 

Fx,Fy» F Z summation of all forces along longitudinal, lateral, and vertical body 
axes, respectively 

F-|,F 2 general terms for longitudinal and lateral perturbation linear and 
angular accelerations (eqs. (15), (16), and (17)) 
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?b 

G 

g 

g n 

H 3 


vector for acceleration of body relative to inertial space, m/sec 2 
(eq. ( 19 )) 

9x6 forcing matrix (eq. (10)) 
scalar magnitude of gravity, m/sec 2 
gravity vector in geographic frame, m/sec 2 
station altitude, m 


I XX> ^YY> -^ZZ 
I XZ 


aircraft moments of inertia about longitudinal, lateral, and ver- 
tical axes, respectively, kg-m 2 


aircraft product of inertia in longitudinal plane, kg-m 2 (JXZ for com- 
puter example in appendix) 

I 1,12,13 inertia products (eqs. (21)) 

aircraft latitude and longitude, rad 

derivatives of aircraft latitude and longitude, rad/sec (eq. ( 24 )) 
latitude and longitude of ground station, rad 


»^s 


Mx,My,Mz rolling, pitching, and yawing moments, respectively 

m aircraft mass, kg 

N random number 

p,q,r aircraft roll rate, pitch rate, and yaw rate, respectively, rad/sec 

Pi >91 > r I measured roll rate, pitch rate, and yaw rate, respectively, rad/sec 


e 

% 

S 


white-noise strength (eq. ( 9 )) 

6 x 1 forcing function vector (eq. (10)) 
range measurement, m (eq. ( 4 )) 
radius of Earth, m 

aircraft altitude relative to mass center of Earth, m (eq. ( 25 )) 
scale factor 


u a» v a> w a measured airspeed components along longitudinal, lateral, and vertical 
body axes, respectively, m/sec 


u a,N» v a,N> w a,N 


nominal airspeed components along longitudinal, lateral, and 
vertical body axes, respectively, m/sec 
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u H» v h» w H aircraft groundspeed components along longitudinal, lateral, and verti- 
cal body axes, respectively, m/sec 

u W> v W’ w W wind velocity components along longitudinal, lateral, and vertical 
body axes, respectively, m/sec 

Ul,vi,wi measured accelerations along longitudinal, lateral, and vertical body 
axes, respectively, m/sec 2 

V n aircraft velocity vector in geographic frame, m/sec 

V n aircraft acceleration vector in geographic frame, m/sec 2 (eq. (18)) 

W wind magnitude, m/sec (eqs. (26) and (27)) 

^a,V,H aircraft airspeed magnitude, m/sec (eq. (26)) 

Xfl aircraft north groundspeed, m/sec 

Xy )H aircraft groundspeed, m/sec (eqs. (12) and (26)) 

x-j,X 2 state perturbations (eq. (15)) 

Y h aircraft east groundspeed, m/sec 

Z H aircraft altitude, m (eq. (25)) 

Z H aircraft vertical velocity, m/sec (eq. (24)) 

AFx,AFy,AF z longitudinal, lateral, and vertical perturbation forces, respec- 
tively in body frame, m/sec 2 (eq. (14)) 

AMy , AMy , AM z rolling, pitching, and yawing perturbation moments, respectively, 

in body frame, rad/sec 2 (eq. (14)) 

Ap,Aq,Ar perturbations in roll rate, pitch rate, and yaw rate, respectively, in 
body frame, rad/sec (eq. (14)) 

ARd range-difference measurement, m (eq. (7)) 

Au a ,Av a ,Aw a longitudinal, lateral, and vertical airspeed perturbation compo- 
nents, respectively, in body frame, m/sec (eq. (13)) 

6l, 6& navigator measured position errors in latitude and longitude, rad 
(eq. (11)) 

<5L,6fi. derivatives of $L and 65, 

(5X perturbation state vector (eq. (11)) 

<5X derivative of perturbation state vector (eq. (10)) 
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6z 

6Z 

e n» e e» e d 

n 

0 az> 0 el 

0Hi4>H»^h 

©h»4 > h»4'h 

6 N><t>N 

yi,u 2 

pj 

Pk 

a 

T 


navigator measured position error in altitude, m (eq. (11)) 
derivative of 6Z 

<5r longitudinal, collective, lateral, and directional control stick 
positions, respectively, m (eq. (14)) 

perturbation control stick vector, m (see fig. 3) 

north, east, and downward platform tilt errors, respectively, rad 
(eq. (11)) 

angle of wind vector measured clockwise from north to tail of vector, 
rad (eqs. (26) and (27)) 

measured azimuth and elevation angles, rad (eqs. (6) and (5)) 

aircraft pitch, roll, and yaw Euler angles, respectively, rad 
derivatives of 0 h> 4>H > and ^H> respectively, rad/sec (eq. (22)) 

measured pitch, roll, and yaw Euler angles, respectively, rad 
nominal pitch and yaw Euler angles, rad (eq. (19)) 
control-stick perturbations, m (eq. (15)) 

measured distance between master station and aircraft, m (eq. (8)) 
measured distance between slave station and aircraft, m (eq. (8)) 
standard deviation 
time constant, sec 


i>o 


initial aircraft groundtrack angle, rad (eqs. (26) and (27)) 


en 


o5 n 

en 

3 b 

ib 


co L 


nb 


u) b 

ib 


skew symmetric matrix for angular velocity of geographic frame rela- 
tive to an Earth fixed frame in geographic-frame coordinates, rad/sec 
(eq. (18)) 

angular velocity vector whose elements comprise matrix S2 n , rad/sec 

en 

angular velocity vector of body frame relative to an inertial frame in 
body-frame coordinates., rad/sec 

angular velocity vector of body frame relative to a geographic frame in 
body-frame coordinates, rad/sec (eq. (23)) 

angular acceleration vector of body frame relative to an inertial frame 
in body-frame coordinates, rad/sec^ (eq. (20)) 


5 



Stability derivative notation : 

be Fy, Fy, F 2 > My, My! or M 2 and A 
v a» w a» P» Q> r > *C. «S. or 5 R 

b,e,i,n body frame, Earth frame, inertial frame, and navigational geographic 
frame, respectively 

H aircraft (specifically helicopter) 

I inertial 

N nominal 

s station 

V vertical heading frame 

W wind 


B * 

3 A 


I trim 
Subscripts : 

a airspeed 


where B can 
can be u. 


“a* 


AVIONICS SIMULATION 
General Flow 

The computer program is composed of many independent modules that represent 
either specific hardware onboard an actual aircraft (such as rate gyros), exter- 
nal navigation aids, or a physical function (such as winds). A simplified block 
diagram of the simulation is shown in figure 1. In general, the helicopter model 



Figure 1.- Simplified block diagram of simulation. 
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calculates and feeds back perfect signals, such as linear and angular veloci- 
ties, accelerations, positions, and angles. The sensor models use as inputs the 
helicopter derived signals and generate appropriate output variables to repre- 
sent the measurements. For perfect measurements, the output is equivalent to 
the input. For models with errors, the output is a corrupted version of the 
input. Each of the three blocks for sensors actually represents several specific 
sensor modules that are described subsequently. 

The outputs from the inertial measuring unit (IMU) sensors feed a navigator 
computer module. The presently available module represents an error model of 
the navigation computer, but the module can easily accommodate actual navigation 
computer equations. 

All of the measurement data feed into an estimator module. A Kalman filter 
is the most likely candidate for the module, although any type of estimator may 
be included. Presently, a dummy model is installed with the output estimate 
being an exact duplicate of the input. 

The guidance computer and external guidance commands are composed of sev- 
eral modules. Guidance commands can be generated internally by stored nominal 
trajectory data and automatic guidance steering laws. External commands are gen- 
erated by the pilot or by air traffic control (ATC). The guidance outputs gen- 
erally represent some type of perturbation error commands and are fed into the 
control computer. The control computer uses these commands to generate equiva- 
lent control stick commands. For the CH-47B helicopter, four control channels 
are used. The model is also capable of accepting input stick commands from the 
pilot . 

The helicopter module consists of several submodules. The main sections 
include actuator and rotor nonlinearities and dynamics, force and moment calcula- 
tions, equations of motion, and an update section to calculate all of the vari- 
ables that are fed to the sensor modules. 


Mainline 

The mainline is essentially an executive that determines the time a module 
is to be executed. Each module is tested to determine the time to update the 
equations, whether it contains differential equations, and if the print option 
has been selected. The iteration frequency for each module is selected by 
either inputting the data or allowing the default value to be chosen. An initial- 
ization subroutine (TIMVARS) uses this frequency data to determine the update 
time increments and the initial starting times. For input frequencies of zero, 
the initial starting times are made very large such that the module is never 
called. The initial starting times for all other modules are set to zero, 
except for the external navigation aids which have selectable start and end 
times . 

In the real-time version, if a module has differential equations, the mod- 
ule is entered and the differential equations are integrated independently of 
the update time. The entire module, including differential equations, is exe- 


cuted when the update time is reached. In the batch program, the differential 
equations are integrated only during the module execution time. 


Inertial Measuring Unit (IMU) Modules 

The IMU modules consist of four major modules , each of which have three mea- 
surements for attitudes, rates, linear accelerations, and angular accelerations. 
The last item is a dummy module with an iteration frequency that presently 
defaults to zero. This module is included to provide the capability of incorpo- 
rating angular acceleration measurements if the need arises. 

The three outputs from the attitude module are the body attitudes for pitch 
6 1 , roll 4* j , and heading ^j, all relative to a local level geographic Earth 
frame. The outputs from the rate gyro and linear accelerometer modules are body- 
axis measurements relative to inertial space. The three rate gyro outputs are 
the pitch rate qj, roll rate pj, and yaw rate r-j-. The three accelerometer 
outputs are the specific force measurements uj, Vj, and Wj in the longitudi- 
nal, lateral, and downward directions, respectively. 

Two types of models are presently offered: a perfect sensor model with the 

output equal to the input and an error model with the output being a corrupted 
version of the input. All three error models have the capability of including 
both bias and random type errors. Both bias and random errors are calculated 
for every sensor by multiplying the output of a white-noise random-number genera- 
tor with Gaussian distribution N by the standard deviation a of the error; 
that is, 

E = Na (1) 

where E represents either type of error. The values of a for both the bias 
and the random errors are input by the user, the default case being zero. Bias 
errors are initialized prior to the operate loop, whereas random errors are 
updated each time the error model is called. 

Scale factor errors are also included in both the attitude-gyro and rate- 
gyro error models. The scale factor S for every gyro is initialized prior to 
the operate loop as 

S = O.OINa + 1 .0 (2) 

Note that all error models simplify to the perfect case if the input a for 
each error is set equal to zero. The capability exists for incorporating sev- 
eral other types of user supplied models. This capability is discussed in the 
section entitled "Preprocessor." 


Other Onboard Sensors 

The four modules included in the other-onboard-sensors block of figure 1 
are air data, groundspeed, heading reference, and absolute altitude. Only per- 
fect sensor models are presently included. The main reason for incorporating 
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these modules is to have the capability of adding models which represent these 
avionics systems with little impact on the overall computer program. The typi- 
cal outputs from each module are defined and are in common storage, but the out- 
puts and their definitions are subject to change, depending upon user prefer- 
ence. The default frequency is zero for each of these modules. 

The presently defined air data outputs are the three body-axis components 
of true airspeed u A d> v A p, and w AD , barometric altitude Z BA r), and angles of 
attack a A D and sideslip 8 AD* Other typical air data outputs that might be 
considered, but are not part of the present module, include indicated airspeed, 
Mach number, and altitude rate. 

The three outputs from the groundspeed module include the north, east, and 
downward components of Earth velocity. A typical onboard sensor that would give 
this type of . measurement is Doppler radar. 

The presently defined heading reference output is the ground-track angle 
relative to north; that is, 


^heading reference = tan 1 -Ji (3) 

X H 

• • 

where Xg and Yg represent the aircraft north and east groundspeeds, respec- 
tively. Another possible heading model that a user might incorporate is for the 
airspeed heading vector. Typical sensors for heading include a magnetic com- 
pass, directional gyro, and gyrocompass. 

The absolute altitude output Z AA is the altitude relative to the ground. 

A radar altimeter and an acoustic altimeter are typical types of sensors. This 
module also has the capability of inputting the height of the terrain directly 
below the aircraft. Presently, the terrain altitude is defaulted to zero, and 
Z AA is the same as Zg A D • 


External Navaids 

Radio communication between ground-based stations and the aircraft provides 
external navigational data. The computer program has the capability of incorpo- 
rating and using simultaneously up to nine different navaid modules in any com- 
bination, plus an additional capability of choosing any one of a possible nine 
models for each navaid. The three navigational measurements presently available 
are for range (DME), elevation and azimuth angles (MLS), and range difference 
(Loran, Omega). Each navigational type has two models from which to select. All 
models have calculations of the actual measurement based upon actual station 
coordinates (latitude, longitude, and altitude), actual helicopter coordinates, 
and noise related errors if the errors are selected. The second model of each 
navaid type also has a calculation of the nominal measurement based upon esti- 
mated station coordinates and desired aircraft coordinates. The nominal measure- 
ment would likely be used in a Kalman filter that is linearized about a nominal 
path. The station coordinates are input data for each of these modules. 
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Distance-measuring equipment (DME) is the range measuring external navaid 
presently installed. A total of 10 stations can be included in any run, with 
any 3 stations operating simultaneously. Call statements to the modules are 
functions of the iteration frequency and the start and stop time for any station 
measurement. Bias and random errors are included in the actual measurement and 
are calculated as described in equation (1). 

The standard deviations for these errors and the time data for each station 
are inputs to the program. The actual calculation range measurement R is 
(ref. 3) 

R = ( R e 2 [(L H - L s ) 2 + (cos 2 L s )0L h -^s> 2 ] + (“ Z H - H s> 2 } + E (4) 

where Ljj, and Z# are the aircraft coordinates (with Z H positive down), 

L s , SL s , and H s are the station coordinates, R e is the Earth radius, and E 
is the noise error. If the aircraft is not within line of sight from the sta- 
tion, a possibility with spherical Earth, or is farther than the maximum DME 
range of 370 km (200 nautical miles), the measurement calculation will not be 
made. The maximum range can be varied by changing a DATA statement. A flag out 
put is available whenever a valid measurement is made. If the print option for 
this module is selected, the program will print any of up to three measurements 
simultaneously, including the time of the last valid measurement for each sta- 
tion. Equation (U) is a good approximation within the DME measurement range. 

For the maximum range of 370 km, the error is approximately 0.15 percent high as 
compared to a more exact formulation that uses spherical trigonometry. 

Microwave landing systems (MLS) are typical angle measuring external 
navaids. The general equation for the actual elevation angle 0 e ^ is 

6 el . tan-’ ~ Z H - "a - 0-5R e r«-H - l s > 2 * (oo» 2 L,H*h - E s ) 2 1 + E (5) 

R e |[(LH - L s ) 2 + (cos 2 L s )(J&h - *3> 2 ]’ /2 

The last term in the numerator is a correction factor since MLS systems gener- 
ally provide angles relative to a coordinate frame that, has the origin at the 
azimuth antenna. The equation for the actual azimuth angle 9 az is 

6 az = tan' 1 (cOS L s )(il H - £s> + E (6) 

e H - L s 

The angles calculated in equations (5) and (6) are applicable to a conical mea- 
suring type system. 

The error E is different in each case. The inputs are similar to those 
for the range measurement, with the addition of a center-line angle for the azi- 
muth measurement. A measurement will not be made if the azimuth angle is out- 
side a 60° boundary from the center line, the elevation angle is not within the 
limits of 1° to 20°, or the maximum range of 55.6 km (30 nautical miles) is 
exceeded. The limits can be changed by modifying a DATA statement. 
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The range-difference models can be used for either Omega or Loran. Inputs 
include position coordinates for both one master and two slave stations 
(ref. 3), iteration frequency, start and end times for a complete measurement, 
and time lags for each of the slave station measurements. The time lags are 
appropriate for Omega where each station has an independent transmission time 
with a constant iteration period of 10 sec. The measurement for each line of 
position is a range difference AR d between the master and a slave as given by 
(ref. 3): 


r AR d = Pj - Pk + E (7) 

where 

Pi = R e cos -1 j^sin Lfl sin L s + cos Ly cos L s cos (5,fj - 5, S )J (i = j,k) (8) 

and Pj is the distance between the master station and the aircraft and is 

the distance between a slave station and the aircraft. A total of two range dif- 
ference measurements are made. Measurement errors include a random error and a 
correlated error. The correlated error is appropriate for Omega and is modeled 
as a first-order time lag to represent a slowly varying bias. Both the time con- 
stant x for the correlated error and the standard deviation a must be input 
to override the default values of zero. The white-noise strength Q is 
(ref. 3) 


q = (9) 

T 

and is calculated in the range-difference initialization module. 


Navigation Computer 

The nine outputs from the navigation computer are the three position coordi- 
nates, three velocity coordinates, and three planform tilt errors when error 
models are included. Presently, three model choices are offered. The default 
model is the perfect case with the output equivalent to the input. The other 
two choices are error models. The second model (first error model choice) incor- 
porates seven states and is representative of most present inertial navigation 
systems. The vertical channel is not included since it is unstable without an 
external reference. The last model includes all nine states and is included to 
provide the capability of future evaluation. In this section, the error model 
description is the nine-state model, since it is the most general case. 

The general error model takes the form of a first-order , nine-state , vector 
differential equation (ref. 4) as follows: 

6$ = F6X + Gq (10) 

where 

= [e n , e e , e d , 6L, 65., 6L, 65,, 6Z, 6z] T (11) 


11 



and 


e n> e e> e d 

6L,6fi,,SZ 

6L,6fl,,SZ 

F 

G 

q 


north, east, and downward platform tilt errors 
latitude, longitude, and altitude rate errors 
latitude, longitude, and altitude position errors 
9x9 system matrix 
9x6 forcing matrix 

6x1 forcing function vector of inertial system errors 


The six inertial system errors consist of the three rate gyro errors and three 
accelerometer errors, all transformed from the body axis to a local level geo- 
graphic frame. The inputs to this module are the nine initial errors which are 
defaulted to zero if the input data are not included. 


Estimator 

The estimator module (see fig. 1) is for a Kalman filter or any other type 
of filter estimator to be installed. Inputs to the estimator are obtained from 
various sensors, external navaids, and the navigation computer. Presently, two 
dummy models are installed. The default model assumes perfect estimates of the 
aircraft position and velocity and of the wind velocity and direction. This 
allows perfect states to be used for guidance and control while nonnavigational 
type errors are evaluated. The second dummy model assumes estimates that are 
identical to the outputs of selected error models, allowing corrupted signals to 
be fed back to the guidance and control system. 


Guidance System 

Guidance commands can be generated in any of four major modules. The first 
module is for storing and calculating the nominal trajectory. This module con- 
tains all of the preflight data and path generating algorithms such as leg 
switching, velocity commands, altitude commands, turning commands, nominal posi- 
tion commands, spiral descent guidance commands , and others. The second module 
is for automatic guidance steering laws. This module contains correction com- 
mands that modify the nominal trajectory as a function of aircraft deviation 
from the commanded path. Typically, this includes modifying the nominal alti- 
tude rate due to deviations in altitude, modifying the nominal bank angle due to 
errors in heading and lateral position, and modifying the nominal acceleration 
or velocity due to errors in velocity and time. The latter is typically a four- 
dimensional time-constrained law. 

The third and fourth sources of guidance commands are the air traffic con- 
trol (ATC) and the pilot. Commands from either of these sources can override 
the earlier guidance commands, the pilot having highest priority. Presently, 
the contents of both these modules are empty. 
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Another function of the guidance computer is to generate error perturbation 
commands for the control system. Presently, only velocity perturbations, based 
upon differences between command velocities and measured velocities, are gener- 
ated. The velocity perturbations are first generated in a geographic coordinate 
frame and then rotated to the body axis for use by the controller. 

Generally, a guidance steering law model is related to a nominal trajectory 
model; for the present discussion, a model pair refers to the combination. 

There are two model pairs presently available. The first pair (default case) is 
for enroute type guidance and can handle data for up to 10 waypoints. Typical 
input data include waypoint position coordinates and desired aircraft velocity 
and altitude. It is also possible to input desired aircraft rate of descent and 
nominal time at each waypoint, although these data are not used in the present 
guidance laws. Separate submodules are used to calculate velocity commands 
(VELCOM), altitude rate commands (CLIMB), and turn commands (LATTURN). A submod- 
ule for calculating the nominal position (PSITION) is also included. Algorithms 
for calculating the correct leg of the trajectory and the times to change vari- 
ous guidance commands are included in submodule LEGSWCH. The only guidance 
steering law in the first model pair is an altitude hold submodule (ALTHOLD) 
that modifies the altitude rate as a function of altitude error. 

The second model pair is for both the enroute phase of flight and for spi- 
ral descent guidance (ref. 5). The enroute guidance is identical to that in the 
first model pair. When the last waypoint is reached, the spiral descent guid- 
ance is automatically switched in. If only one enroute point is included in the 
simulation run, the program will start with the spiral descent guidance laws. 

The input data required are the coordinates of the landing pad and spiral cen- 
ter, airspeed, spiral turning radius, and turn direction. The last three items 
have default values. 


Control Computer 

The control-computer module is used for any of the research control laws to 
be tested. Typical inputs include attitude and rate measurements, guidance com- 
mands, and velocity error perturbations. The outputs are four stick commands 
for the helicopter control inputs. The default model presently installed has a 
velocity control law. 


Winds 

Three wind models are presently offered. The default model is the zero 
wind case and sets all wind variables to zero. The second model is for either a 
constant horizontal wind vector or a wind shear. The wind shear allows a varia- 
tion of the horizontal wind vector with altitude. Up to four data points can be 
included for both wind velocity and wind direction relative to north. Straight 
line interpolation is used, with altitude the independent variable. If only one 
data point is entered, a steady wind will be used, and if no data points are 
entered, a zero wind will be calculated. The vertical wind is assumed zero in 
this model. 
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The third wind model contains the contents of the second model plus vari- 
able winds in both horizontal and vertical channels. The variable winds are 
modeled as exponentially correlated winds with specified correlation distances 
dty and standard deviations a, plus random gusts (ref. 3). The white-noise 
strength Q is calculated as shown in equation (9), where a separate g is 
input for the horizontal velocity, horizontal direction, and vertical velocity. 
The three correlation distances have defaulted values but can be overridden by 
input data. The effective time constants t for the exponentially correlated 
winds are a function of d^ and the groundspeed Xy , H 5 that is, 


T 



( 12 ) 


The random gusts for both horizontal velocity and vertical velocity are calcu- 
lated as in equation (1). The third model reduces to the second model by set- 
ting the input standard deviations to zero, and both models reduce to the 
default model with all input data defaulted to zero. The iteration frequency of 
the wind module is automatically set to the aircraft model iteration frequency 
and cannot be adjusted. 


AIRCRAFT SIMULATION MODEL 

A Boeing Vertol CH-47B helicopter is being used as a tool to evaluate 
advanced avionics research concepts as part of the automatic avionics program 
(ref. 1). A model of the CH-47B is included in this computer program to evalu- 
ate the avionics concepts prior to flight. One model is presently offered, 
although the preprocessor capability and modular concept allow other models to 
be added in the future. A linear perturbation model, for computing aerodynamic 
forces and moments, has been selected over a nonlinear model (ref. 6) to mini- 
mize computation time. A unique technique for updating the trims and deriva- 
tives, to maintain small perturbations over the complete flight regime, is 
described in this section. The kinematic equations of motion are fully coupled 
nonlinear equations. 

The helicopter model consists of several submodules as shown in figure 2. 
Typically, control inputs originate from the control computer, including a sta- 
bility augmentation system, and from the pilot. These control forces are inputs 
to the actuator and rotor section which includes nonlinearities and dynamics. 

The control-stick trim nominals are combined with the equivalent total control- 
stick inputs to generate perturbation stick commands. The airspeed submodule 
combines airspeed trim nominals, winds, and actual helicopter velocities to cal- 
culate perturbation airspeeds. Perturbation airspeeds, perturbation stick com- 
mands, and perturbation angular velocities are combined with the stability and 
control derivatives to calculate the total perturbation forces and moments act- 
ing on the helicopter. The total helicopter forces and moments are calculated 
by adding nominals and perturbation forces and moments. The equations of motion 
contain 12 differential equations for linear and angular accelerations and veloc- 
ities, and take into account spherical Earth. In the update module, all of the 
helicopter variables are updated for the next time period. Each of the helicop- 
ter submodules is described in more detail in the following sections. 
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Figure 2.- Block diagram of helicopter module. 


Actuator and Rotor Section 

An analysis has been made of the CH-47B control linkage system, from stick 
inputs to rotor motion , and the dominant system time constants and nonlineari- 
ties are taken into account. (See fig. 3.) One change in the research CH-47B 
model from a standard model is the addition of electrohydraulic actuators 
between the control sticks and the transfer bell cranks. The actuator position 
limit (+5.08 cm from ref. 6) is the limiting mechanical displacement between the 
actuator input and rotor blades. It is possible for the pilot's control stick 
to be physically limited before the actuator limit occurs, but the stick limits 
would be modeled in the pilot module rather than the helicopter module. The 
dominant-system rate-limit constraint comes from the electrohydraulic actuator 
(+7.62 cm/sec from ref. 6). Neither the upper boost nor the lower boost actua- 
tors (ref. 7) will reach their respective rate limits if the electrohydraulic 
actuator rate limit is not reached. The time constants of both the electro- 
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stick command, 


Figure 3.- Block diagram of CH-47B actuator and rotor section. 

hydraulic actuators (0.0125 sec from unpublished data) and the lower boost actua- 
tors (0.0024 sec from ref. 8) are extremely small and are ignored. The upper 
boost actuators have an approximate second-order response with a damping ratio 
of 0.55 and a natural frequency of 50 rad/sec. (See ref. 6.) Their effect on 
the total system dynamics is minor when compared with the rotor flapping dynam- 
ics which has a damping ratio of approximately 0.6 and a natural frequency of 
24 rad/sec (from unpublished data). The blocks in figure 3 labeled "Gain" and 
"1/Gain" are included to represent actual physical quantities. The gains for all 
four channels are selected such that a full electrohydraulic actuator motion 
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corresponds to a full pilot control stick throw. The second gain block is 
included to maintain the same total system steady-state gain used in calculating 
the tr ims and derivatives (ref. 2) . The output from the rotor dynamics repre- 
sents an equivalent total stick position. This value is summed with the corre- 
sponding nominal control stick position to get the perturbation stick command 

3 t . 
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The three body-axis airspeed perturbations 
lated in the airspeeds submodule as 
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(13) 


where ujj, v^, and wy are the actual aircraft body-axis velocity components, 
U W> V W> and w W are the three wind components in body axis, and u a fj, v a,N» 
and w a ^ are the nominal airspeed trim components in body axis. The approach 
for calculating the trims is described later in this section. 


Forces and Moments 

Perturbation linear accelerations AF^/m, APy/m, and AF z /m and angular 
accelerations AMy/Iyy, AMy/Iyy, and Am z /I zz are calculated in the forces and 
moments submodule by linearly combining the aerodynamic and control force and 
moment contributions as follows: 
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The aerodynamic driving terms are perturbations in airspeeds Au a , Av a , 
and Aw a and angular velocities Ap, Aq, and Ar, and the control force inputs 
are perturbations in the four control stick inputs 8 b > 8c > Ss> and Sr. Equa- 
tion (14) is partitioned to separate the longitudinal and lateral equations 
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The A matrix terms A-| , A 2 , A3, and A4 represent the stability deriva- 
tives, the B matrix terms B-] , B2, B3, and B4 represent the control deriva- 

tives, the state perturbations are represented by x-| and x 2 > and the control 
stick perturbations are represented by p-j and p 2 . Equation ( 15 ) is shown in 
the most general form with coupled aerodynamics. Generally, the coupling deriva- 
tives A2, A3, B2, and B3 have negligible effect relative to the uncoupled 

derivatives A -j , A4, B-j , and B4. Only the uncoupled derivatives are being 

used in the present model to save storage and computing time. All of the deriva- 
tives are available (ref. 2 ) and can be incorporated. Derivatives for the com- 
plete flight regime are stored in a stability derivative module (STABDER), and 
the technique for updating the derivatives shown in equation ( 14 ) will be 
described. The derivatives are updated at regular intervals to maintain small 
perturbations from nominal trim conditions. One second appears to give suffi- 
ciently good results. 


Kinematic Equations of Motion 

This section describes the kinematic-equations-of-motion submodule; 12 dif- 
ferential equations are included to describe the total forces and moments acting 
on the vehicle, Euler rates, and Earth referenced velocities. The equations are 
for a spherical, nonrotating Earth. The linear acceleration and velocity equa- 
tions are referenced to a local level geographic coordinate frame whereas the 
angular acceleration and rate equations are referenced to the body axis. 

The vehicle acceleration vector fa relative to the local level geographic 
frame is 



- OP fa 

en 


( 18 ) 


where C5 is a matrix transformation from body-axis frame to geographic frame, 
fa is the gravity vector in the geographic frame, ftn is a skew symmetric 

matrix for the angular velocity of the geographic frame relative to the Earth 
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fixed frame, V n is the aircraft velocity vector in the geographic frame, and 
fb is the specific force vector in the body-axis frame (ref. 9). The specific 
force vector is calculated by summing the perturbation linear accelerations 
AF x /m, AFy/ai, and AF z /m generated in the force and moment submodule with the 
gravity dependent nominal trim terms; that is, 
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(19) 


where % and are the trim pitch and bank angles that are generated in sub- 

module STABDER . The components of f b are fed back to the linear accelerometer 
module discussed in the IMU section . 


The body-axis acceleration vector 
approximated as 
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relative to inertial space is 
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The perturbation angular accelerations AM x /I xx , AMy/IyY* and ^Z^ZZ calcu- 
lated in the forces and moments submodule are proportional to the total moments 
acting on the helicopter (ref. 2). The inertia products I-| , l 2 > and I 3 are 
functions of the moments of inertia I XX , Iyy > and ^ZZ and product of inertia 
I xz in the helicopter XZ-plane as follows: 


I-, r. I XX I ZZ 

I xx I zz - dxz ) 2 


( 21 a) 


I = I 1 I XZ (21b) 

I ZZ 

lo = I 1 I XZ (21c) 

5 J xx 

The three integral components of equation (20) are fed back to the rate gyro mod- 
ule discussed in the IMU section. 
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The Euler angle rates , using a yaw 
order are 
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The three integral components of equation (22) represent the Euler angles that 
relate the body frame to the local level geographic frame and are fed back to 
the attitude module discussed in the IMU section. The angular velocity vector 
->b 

tonb of the helicopter relative to the local level geographic frame is related 
to the difference between the body rates relative to inertial space and the angu- 
lar velocity of the local level frame relative to the Earth fixed frame; that 
is, 
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(23) 


The rate of change of the helicopter latitude 
tude Zjj relative to the Earth is 
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(24) 


where Xj^ is the north component of groundspeed, is the east component of 

groundspeed, and (with Ry its derivative) is the vehicle altitude rela- 

tive to the mass center of the Earth. The altitude R H is a function of Earth 
radius R e and helicopter altitude Z H above the Earth surface as 


R H = R e " Z H 


(25) 


The Earth radius is assumed to be a constant of 6378.1631 km (20 925 732 ft). 

The integral of equation (24) gives the Earth related coordinates. In this pro- 
gram, latitudes are positive in the northern hemisphere and longitudes are nega- 
tive in the western hemisphere. 


Updates 

In this section, all of the helicopter variables are updated for the suc- 
ceeding time period and for distribution to all of the other modules in the sim- 
ulation. Body referenced helicopter variables include both linear and angular 
accelerations, linear and angular velocities, attitudes, and angles of attack 
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and sideslip. Earth referenced helicopter variables include linear accelera- 
tions, velocities, and positions. Coefficients of a direction cosine matrix are 
also available for use in other modules. 


Helicopter Initialization 

The input data for the helicopter module consist of the iteration frequency 
(default value is 32 iterations per second) and three sets of three initial 
errors (default values are zero). The first set is referenced to the first way- 
point location and includes along- track position error, cross-track position 
error, and altitude error; the second set is referenced to the desired veloci- 
ties and heading at the first waypoint and. includes groundspeed error, ground- 
track error, and altitude rate error; and the third set is body referenced and 
includes pitch attitude error, bank attitude error, and heading error. 


The airspeed vector is calculated by summing the groundspeed and wind vec- 
tors. The airspeed magnitude x a ,V,H is 


where 


^a,V,H = k [(^V,H) 2 + W 2 + 2Xy >H W cos ^ “ ^o^] 
+1, (Xv^ + W cos (ri - ^ 0 ) > 0) 
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and the heading of the airspeed vector \pfj relative to north is 
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(27) 


where Xy H is the groundspeed, W is the wind magnitude, \|j 0 is the initial 
ground track, and n is the wind direction relative to north. Any initial head- 
ing error is added to equation (27). 


The airspeed X a y ^ calculated in equation (26) is used in the stability 
derivative and trim lookup table to find the equilibrium pitch and bank atti- 
tudes. Initial pitch and bank attitude errors are then added to the equilibrium 
values to determine the initial helicopter attitudes. A complete set of deriva- 
tives and trims is also determined at this time. Other calculations made in 
this module include the three specific forces and the three body rates relative 
to inertial space. Finally, the 12 differential equations described in the 
equations-of-motion section are initialized. These equations include the three 
groundspeed components, three body rates, three attitudes, and three position 
components . 
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Stability Derivatives (STABDER) 

The purpose of STABDER is to determine the stability and control deriva- 
tives and trim nominals during flight, at any operating point in the flight 
regime. The module is essentially composed of two parts; the first part deter- 
mines the operating points at which to calculate the derivatives and trims, and 
the second part is the calculation of the derivatives and trims. The flight- 
regime of the CH-47B helicopter includes airspeeds between -40 and 160 knots and 
rates of descent between +10.16 m/sec (+2000 ft/min). The nominal operating 
point is a function of both the commanded helicopter velocities (airspeed and 
rate of descent) and actual helicopter velocities and is always maintained 
within predefined limits of the actual helicopter velocities. When the differ- 
ence between the actual and commanded velocities is less than the predefined 
limits, the commanded velocities are chosen for the nominal operating point. 

The assumption is that for a stable vehicle the helicopter velocity will con- 
verge to the commanded value. When the difference between the actual and com- 
manded velocities is greater than the predefined limits, the nominal operating 
point is chosen as the sum of the actual velocity and the predefined limit. The 
airspeed limit is chosen as 9-6 knots and the rate of descent limit is 1.5 m/sec 
(300 ft/min). These limits are defined in DATA statements and can easily be 
changed. The operating point is determined every second, and is also set in a 
DATA statement. If the nominal operating point does not change during the 1-sec 
interval, the previously calculated derivatives and trims are used. 

For the uncoupled system, there are 48 derivatives and trims to calculate 
each time the nominal operating point is changed. Furthermore, for each deriva- 
tive and trim there are 99 data points; this includes 11 airspeeds at 9 rates of 
descent. In order to save storage, third-order curves are fitted through each 
of the 11 airspeed data points. This approach allows the storage to be reduced 
from 4752 -jo to 1 728 -| q locations. For each derivative and trim, the third-order 
algebraic equation is first calculated as a function of airspeed; linear inter- 
polation is then applied, using rate of descent as the independent variable. 
Initial condition and step input responses of the open-loop helicopter model 
show that the third-order curves have approximately the same dynamic responses 
as those obtained by using all of the derivatives . 


PREPROCESSOR 

The simulation program has the capability of selecting various combinations 
of models for any run from a choice of many different model types, thereby mini- 
mizing the in-core storage and computation time. The feature that allows this 
capability is the preprocessor. In general, the preprocessor is a FORTRAN pro- 
gram that performs modifications on the main program prior to compilation and 
execution of the main program. The most obvious function of the preprocessor is 
to decide upon the appropriate models to be used and to write the correct CALL 
statements into the computer program. Model selection is accomplished by a com- 
bination of defaults and input data cards. Specific functions include: (1) 

writing subroutine calls within the mainline and within other subroutines, 

(2) choosing appropriate initialization models and writing the correct CALL 
statements within the initialization routine, (3) writing the necessary COMMON 
blocks for on-line plotting, (4) calculating the total number of differential 
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equations to be solved in the simulation and writing the COMMON blocks with 
proper array dimensions, (5) writing specific CALL statements to the real-time 
system for pseudo real-time runs, and (6) printing an alphabetical table of all 
subroutines used in a specific run and the total number of differential equa- 
tions in the simulation. 

The means for determining the correct subroutine CALL statement is 
through a module integer number (KEY 1) and a model integer number (KEY 2). 
Every major selectable module and submodule in the simulation program have an 
assigned KEY 1 integer (10, 20, 30, etc.). The model choices within a module 
have assigned one-digit numbers (0, . . ., 9), allowing a possibility of up to 
10 model choices. The sum of KEY 1 and KEY 2 relates to an array SNAM and 
defines the particular subroutine to be called. 

A general flow chart is illustrated in figure 4. Default values are first 
stored in all arrays and then are overridden by user input data cards. Input 



Figure 4.- General flow chart for preprocessor. 


23 









data cards are only required if optional subroutines are to be chosen. The 
unprocessed main program has the word "FLAG" in columns 1 to 4, wherever a 
selectable subroutine is to be called. The preprocessor checks each statement 
in the main program. If the word "FLAG" or one of the other special indicators 
(not defined here) for function 3, 4, or 5 defined previously is not encoun- 
tered, the statement is written onto the file to be executed. If the word 
"FLAG" is encountered, the preprocessor checks the very next statement for an 
index integer number. The index number corresponds to the correct module (or 
submodule) to be called. Using this index number along with KEY 1 and KEY 2 
data, the SNAM array location, corresponding to the correct subroutine to be 
called, is calculated. The preprocessor then writes this subroutine call onto 
the file to be executed. This process is continued until an end of file is 
detected. 

Input data cards are only used when the default cases are not desired. 

Each data card includes three numbers: the first corresponds to the KEY 1 inte- 

ger corresponding to the module being selected, the second corresponds to the 
KEY 2 integer corresponding to the model being selected, and the third integer 
is the number of differential equations in the model. For example, a data card 
with integers 200, 1, and 7 indicates a local level navigator module, model num- 
ber 2, with 7 differential equations. Two problem examples for the preprocessor 
are illustrated in the next section of this paper. Presently, 39 modules have 
been defined for the batch program and 44 modules, for the real-time version. 

The preprocessor has the capability of writing subroutine calls within 
subroutines as well as the mainline. The process is the same as previously 
described. However, there are additional considerations involved. For example, 
a subroutine called within another subroutine need not be added to the subrou- 
tine library if it has already been called in the mainline or in a previous sub- 
routine. The preprocessor insures that a subroutine will only be added to the 
final library once, regardless of the number of times it has been called. The 
preprocessor will handle a subroutine call depth of up to four levels. A tree 
structure is used to keep track of the call depth levels and to determine which 
subroutines have been called and added to the library. 


PROGRAM ORGANIZATION 
Deck Structure 

The deck configuration consists of five records plus the system control 
cards, as shown in figure 5. The first three records are used for modifications 
to key parts of the overall program; the first record is for modifications to 
the stability derivative curve coefficients, the second record is to modify the 
mainline and subroutine library, and the third record is to modify the preproces- 
sor. The preprocessor input data cards, discussed in the preprocessor section 
and illustrated in the problem examples section, are located in the fourth 
record. The fifth record is for the main program input data cards, which are 
illustrated in the problem examples section. 
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System control cards 


EOR 

Mods for stability derivative curve coefficients 


EOR 

Mods for mainline and subroutine library 


EOR 

Mods for preprocessor 


EOR 

Preprocessor data cards 


EOR 

Main program data cards 


EOF 

Figure 5.- Control card deck structure. (EOR is end 
of record and EOF is end of file.) 


Problem Examples 

This section contains some examples illustrating the input data cards for 
the preprocessor and the main program. The preprocessor data are field oriented, 
whereas all main-program input data are in the form of NAMELIST statements. The 
only requirement is that all NAMELIST cards be present in the data deck, even 
though all of the FORTRAN variables defined under a NAMELIST statement have 
defaulted values. Trajectory type data are the only input data that must be 
submitted. 

The first preprocessor example illustrates the default option and shows the 
minimum input data required. The example case is to use all default modules. 

The preprocessor data deck consists of one card with a ”0” in column 5, indicat- 
ing that no other cards are required. 

The second preprocessor example represents a more typical case. It is 
desired to have error models for the navigator (LLNAV2), rate gyros (RATEG2), 
and linear accelerometers (LINAC2) and to use the second wind model (EXTWIN2). 
The subroutine call names are in parentheses. All other modules are defaulted 
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to the first model case. The data deck for this example is illustrated in / 
figure 6. The first card contains a "4" in column 5, indicating four data/cards 
are to follow. After the first card, the names of each model can be included in 
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Figure 6.- Data deck for second preprocessor ekample. 

\ 

columns 1 to 10 for user identification, since the preprocessor does not recog- 
nize these columns. The three numbers on each card correspond to the module num- 
ber (columns 11 to 15), the model number selected (columns 1 6 to 20), and the 
number of differential Equations (columns 21 to 25). All numbers are right jus- 
tified. Alphabetical ,brder is the most efficient configuration, although any 
order can be used. This example will have 19 differential equations, 7 for the 
navigator and 12 defaulted differential equations for the helicopter model. 
Presently, there are 4-4 module numbers reserved, including 7 for the real-time 
program. / \ 

The main /program example uses the four modules described iri the second 
preprocessor /example. In particular, a trajectory is defined, initial errors 
are defined for the helicopter and the various error models, wind\data are 
described ,/ print frequencies are defined, and print data are defiried for several 
modules. The specific example is as follows: \ 


to the first model case. The data deck for this example is illustrated in 
figure 6. The first card contains a "4" in column 5, indicating four data cards 
are to follow. After the first card, the names of each model can be included in 


RATEGE 

£40 

1 

0 

'LLNAV2 

£00 

1 

7 
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190 

1 

0 
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80 

1 
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' 4 
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Figure 6.- Data deck for second preprocessor example. 

columns 1 to 10 for user identification, since the preprocessor does not recog- 
nize these columns . The three numbers on each card correspond to the module num- 
ber (columns 11 to 15), the model number selected (columns 16 to 20), and the 
number of differential equations (columns 21 to 25). All numbers are right jus- 
tified. Alphabetical order is the most efficient configuration, although any 
order can be used. This example will have 19 differential equations, 7 for the 
navigator and 12 defaulted differential equations for the helicopter model. 
Presently, there are 44 module numbers reserved, including 7 for the real-time 
program. 

The main program example uses the four modules described in the second 
preprocessor example. In particular, a trajectory is defined, initial errors 
are defined for the helicopter and the various error models, wind data are 
described, print frequencies are defined, and print data are defined for several 
modules. The specific example is as follows: 


1. Define the nominal trajectory: 


Waypoint 

Latitude , 
deg 

Longitude , 
deg 

Altitude, 
m (ft) 

Velocity, 
m/sec (ft/sec) 

1 

40.64 

-74.067 

-426.7 (-1400) 

30.48 

(100) 

2 

40.64 

-74.04 

-335.3 (-1100) 

24.38 

(80) 

3 

40.66 

-74.025 

-243-8 (-800) 

19.8 

(65) 


2. Assume initial helicopter errors: 

Along track error, m (ft) -30.48 (-100) 

Altitude error, m (ft) 3.040 (10) 

Velocity error, m/sec (ft/sec) 0.9144 (3) 

3. Assume a linear accelerometer bias error of 0.003048 m/sec^ 

(0.01 ft/sec^). 

4. Assume a rate gyro scale factor error of 2 percent. 

5. Assume the navigator has a north platform misalinement angle of 1 yrad. 

6. Assume a wind shear with two data points: 


Point Altitude, Wind direction, Wind velocity, 

m (ft) deg m/sec (ft/sec) 

1 0 ( 0 ) 180 0 ( 0 ) 

2 304.8 (1000) 235 6.1 (20) 

7. Assume all default frequencies except: 

Accelerometer module, per sec 32 

Navigator module, per sec 32 

Printout, per sec 0.1 

Maximum run time, sec 20 


Obviously, the run time for this example will complete only a very minor 
portion of the desired trajectory. 

8. Print the following data: 

All initialization data 
Rate gyro outputs 
Linear accelerometer outputs 
Navigator outputs (option 2) 

Guidance commands (option 1) 

Controller outputs 
All helicopter outputs 

Figure 7 illustrates the input data cards for this example. Eighteen 
NAMELIST data cards are shown; the real-time version has three additional input 
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Figure 7.- Data deck for main program example. 
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Figure 7.- Data deck for main program example. 
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cards. The NAMELIST card DECNTS describes the number of differential equations 
for each module in the simulation run, and the NAMELIST card PVAR describes the 
print variable options for each module. All other data in figure 7 can be corre- 
lated with numbers in the example. The U.S. Customary Units are used in the com- 
puter program. A computer printout for the second preprocessor example and the 
main program example is shown in the appendix. 

The computer example was run on a Control Data Corporation (CDC) Cyter 175 
computer, with a CDC FORTRAN FTN compiler (version 4.6), Version 1.1 of the NQS 
operating system was in use. The actual central processor times for compiling, 
loading, and executing the various sections of the computer program are summa- 
rized as follows: 

, , . 2.043 
. . . 2.951 
. . . 1,086 
, . • 7.959 
. . . 3-819 

The actual storage needed for this example is 50 600 octal words. 


Compile preprocessor, sec ...... 

Load and execute preprocessor, sec 
Compile main program, sec ...... 

Compile subroutine library, sec . . . 
Load and execute complete program, sec 


CONCLUDING REMARKS 

A user-oriented computer-simulation program has been developed to allow for 
the evaluation of individual avionics systems in conjunction with a VTOL air- 
craft model and other types of VTOL avionics, systems. After a model has been 
developed by any user, the subroutine is incorporated into the system library 
for use by all of the researchers. The input data for the main program are 
defined in NAMELIST statements. A series of defaults has been incorporated to 
minimize the amount of input data needed for any simulation run. 

The modular approach allows many different models to be developed and 
installed into the subroutine library, with minimum impact on the total system. 
Individual modules for all of the basic avionics sensors allow a large choice of 
measurements and flexibility. 

A preprocessor has been developed to allow a large selection of models. 

For any module, a possibility of up to 10 models exists. Default models allow a 
minimum of input data cards for the preprocessor. 

A real-time program is maintained as nearly similar as possible with the 
batch program. By using the preprocessor capability, a dummy model with zero 
input frequency can be selected by input data cards; thereby, the software is 
eliminated from storage. The software can be replaced by actual hardware that 
operates through the electrical input/output channels. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton , VA 23665 
February 9, 1977 
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APPENDIX 


COMPUTER PRINTOUT FOR PROBLEM EXAMPLE 1 ' 

This appendix contains the computer printout for the problem example. The first paige of the 
printout lists all of the subroutines selected by the preprocessor and indicates the total number of . 
differential equations for the example. All other pages are related to the main program. Initializa- 
tion data are included on the first 3-1/2 pages; this includes the iteration frequencies of all modules, 
trajectory data, helicopter data, trim data, and all error model initialization data,. The completion of 
the printout contains data from the modules selected at the iteration frequency. 


VALT/FID BATCH SIMULATION 


THE FOLLOWING SUBROUTINE' U" c U'Fn in THIS RUN: 


ABSALT 

ACTOYN 

ACTPHT 

A T°n AT A 

AIRSPD 

CHSIIR 

CLIMB 

COLWTM^ 

L 

ETRAN 

GSPEEO 

GUI DPR T 

GUPDATC 

HCfjD CP 

HE L CON 

IATCM01 

ICnNTRl 

IPXWTM? 

inrmio i 

IGIDPT1 

ILINAC2 

ILLNAV2 

I NTE p 

THTO A J 1 

IPILnTl 

TXNAVA 

IXNAV5 

IXNAVf' 

TYNAW7 

IXNAV8 

LLNAV2 

M E ANW I N 

N AVFO’n 

nhtpaj i 

NR AN2 

POSFST 

POSLIM 

psitthn 

o ATrr,? 

PATELIM 

5PGUI0 

SPTRALl 

SPTRGTn 

q o u T N n 

STABDER 

X N A V 1 

XNAV2 

XNA V** 

YM A V 4 

XNAV5 


ZSTAB 


ALTHOLO 

AMOVE 

ANGAC 

ATCMOD 

BODYEOS 

EXTWIN2 

FANDM 

GCOUPIE 

GENEOS 

GSLAW1 

HELICOP 

HYSTER 

IARSAL1 

IAIR0A1 . 

I ANGAC 1 

IGSLAW1 

IGSPED1 

IHDREF1 

IHELICP 

IKALMN1 

IRATEG2 

IRIGY01 

IXNAV1 

IXNAV? 

IXNAV3 

IXNAV9 

kalmani 

LATTURN 

tEGSWCH 

LIN AC 2 

ONPLOT 

OUTPUT 

PILOT 

BLOTITl 

RLTSAV1 

RDINPT 

R IGYRD1 

ROTDYN 

SPENTRY 

SPEXIT 

SVARDES 

TIMVARS 

TRKDTST 

UPDATE 

VELCOM 

XNAV6 

XNAV7 

XNAVB 

XNAV9 

ZHELICP 


THE TOTAL NUMBER OF DIFFERENTIA) 'OUATIDNS IN THIS RUN IS: 19 




INITIALIZATION SECTION 


HELICOPTER INERTIAS ( PL UG-F fct fc.) 


IX* 

•3716300E+05 


irr 

•2017500E+06 


II? 

• 1 QQ05Q0£+0fi 


FREQUENCIES IN 


ITERATIONS P co « p niA 


HELICOPTER 

RATE INTEGRATING GYRO 
RATE GYRO 

LINEAR ACCELEROMETER 
ANGULAR ACCELEROMETER 
AIR DATA 
GROUNDSPEED 
HEADING REFERENCE 
ABSOLUTE ALTTTUOE 
EXTERNAL NAV. AIO 1 . 
EXTERNAL NAV. AID 2 
EXTERNAL NAV. AID 3 
EXTERNAL NAV. AID A 
EXTERNAL NAV. AIO 5 
EXTERNAL NAV. AID 6 
EXTERNAL NAV. AID 7 
EXTERNAL NAV. AID B 
EXTERNAL NAV. AID Q 
LOCAL LEVEL NAV. 

KALNAN FILTER 
NOMINAL TRAJECTORY 
GUIDANCE STEEPING LAW 
AIR TRAFFIC CONTROL 
PILOT TRAJECTORY 
GUIDANCE PERTURBATIONS 
GUIDANCE COUPLER 
CONTROLLER 
PRINTOUT 


• 3 7O0DF +02 
.12O00B+02 
.3PM00F+02 
.i?000 c +02 

D. 

0. 

o. 

n. 

0. 

0. 

n. 

o. 
o. 
o. 
o . 
o. 
o . 

0. 

.3P0DDF+02 

.iKnnoR+02 
.poooor +oi 

.PDOOOE «01 
0. 

,’ooooc +oi 

.1 BOOOF +02 

• I GOOQF +02. 
.1BO0DF+0Z 

.1 nonoF +oo 


MAXIMUM RUN TIMF = .POnnOF+O? +FCONOS 

INTEGRATION TIME INTERVAL IP . 31 2 50E-0 1 SECONDS 

RANOOM NUMBER GENERATOR SEFD = Ba.OHO 


JX2 

. 1A63200E+05 
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HELICOPTER DATA 


geographic FRAME variable 


LATITUDE (DEG) 
.A0ft4000E+02 

LPNO y t»(of { np Q | 
-,7406736' : +02 

ALTITUDE (FT) 
-. 1390000E+04 



N. VEL. (FT/SEC) 
-•1766578E-12 

E • Vt=t . f FT/Scr ) 

• 1 030000F +03 

TOTAL VEL • (FT /SEC 1 
•1030000E+03 

GROUND TRACK { 0 EG ) 
• 90Q0000E+ 02 

DESCENT FATE (FT/SEC) 

0. 


BODY AXIS VARIABLES 


LONG. VEL. ( F T/SFC ) 

1 ATFRAL '' c L . (FT/SEC ) 

VERT. VEL. 

(FT/SEC 1 

• 101899RE + 03 

-.1 3560ft7c + o? 

.6445725E+01 


PITCH ANGLE IDEG) 

RANK ANGLF (DEG) 

YAW ANGLE 

(DEG) 

• 36D2560E+01 

-.32*570*1^+00 

. 97544 3 3E + 02 


ANGLE OF ATTACK (DEG) * 

, 366^801 c +01 



ANGLE OF SIDESLIP (DEG) 

* -.7QR4R7 7E-01 




TRIM DATA 


groundspefds airspeeds 

U NOM (FT/SEC) V NOMtFT/^n V NOM(FT/SEC) UA NOM(FT/SEO VA NOM(FT/SEC) WA NOMIFT/SEC) 

.989? 2336+02 - . 1 3560? 7F + 0? .*416447E+01 .8421890E+02 3 1 74507E-01 .5528B36E*01 

ATTITUDES 

THETA NOM (DEG) PHI NOM <0^1 
• 375604P F +0 1 - . 3 20 973 ?F +00 


RATE GYROS 


D 


ijr u 


ROLL 


YAW 


P I A S (RAD. /SEC.) 
RANDOM (RAD. /SEC.) 
SCALE FACTOR 


0. 

0. 

, I 0003 9 A P ♦ 0 1 


0. 

0. 

. 1018096E+G1 


0. 

0. 

• 9507966E+OO 


UJ 

U) 


r 


u> 

x= 


LINEAR ACCELEROMETERS 


LONC.. LATERAL VERTICAL 

BIAS (FT. /SEC. SO.) - . U AROA7F-01 -. I213770E-01 .21722R5E-02 

random (ft./sfc.so. ) n. o. o. 


LOCAL LEVEL NAVISATOR 


MISALIGNMENT 

ERRORS 


POSITION 

ERRORS 


VELOCITY 

ERRORS 


north 

EAST 

nnvN 

CAT 

LONG 

ACT 

LAT RATE 

LONG RATE 

VERTICAL 

(RAO) 

(RAO) 

< PAO) 

(RADI 

(RAO) 

(FT) 

(RAD/SEC! 

1 (RAD/SEC) 

(FT/SEC) 

.10000E-05 

0. 

o. 

0. 

0. 

0. 

0. 

0. 

0. 


RATE GYRO OUTPUTS 

PI • -.37517E-06 

01 * -.A9A28E-05 

RT • -.AI083E-05 

TIMF • 0. 


LINEAR ACCELEROMETER OUTPUT* 

UDOTI ■ . 20977E+01 
VDOTI ■ • 1 7235E + 00 
WOOTI » -.3212BE+02 


TIME 


0 
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LOCAL LEVEL NAVIGATOR 


MISALIGNMENT 

ERRORS 


POSITION 

FPRORS 


VELOCITY ERRORS 


NORTH 

EAST 

Ohwn 

LAT 

LONG 

ALT 

LAT RATE LONG RATE 

VERTICAL 

(RAD) 

(RAO) 

f 

(RAD) 

(RAO) 

(FT) 

(RAD/SEC) (RAD/SEC) 

1 FT/SEC ( 

• 10000E-05 

0 . 

0 , 

0 . 

0 . 


I 0 . 0 . 

I 


FINAL GUIDANCE COMMANDS 
TIME ■ 0. SEC 

XOOTGSC (FT/SEC) * .qoMOOF+n’ PSIHPC (DEG) » .90000E+02 7D0TGSC (FT/SEC) ■ .B6Z50E+00 

■ZGSC (FT),* -.13997E+DA 

. RHIGSC (DEG l * 0. °S I GSC (OEGI * .90000E+OZ 


CONTROLLER OUTPUTS (INCH^R) 



OCC(l) 

DCC l 2 ) 

nrr f 3 i 

DCC (A) 

.5B3Z7E+00 

- . 23786F+00 

• -10400E + 00 ■ 

22A99E +01 


TIME *" 0. 



U> 

(J1 
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If 


CO 

O' 


HELICOPTER VARIABLES AT TT“F ■ 0. 

POSITIONS 


XH 

LATH 

■ 0 • 

- .406<iOOF*0? 

YH 

IHNH 

S 

a 

0. 

-•740674E+02 

NAV. GROUNOSP6EOS (FTASEFI 




XOOTH 

* -.176658F-12 

vonjH 

a 

•103000E+03 

BODY VELOCITIES (FT/SECI 




UH 

* •101900F<Mtt 

VH 


-.135601E+D2 

LINEAR ACCELERATIONS (FT/SFF.S9'. 




OUDOTH 

* 0 . 

vrv«TH 


Q« 

EULER ANGLES (DEG ) 




THEM 

-■ a 366 256F +01 

PHTH- 

a 

-.325708E+Q0 

BODY RATES 

(RADI ANS/S EC )' 




PH 

= — a 375029F— Oft 

OH 

a 

-.A85996E-05 

ANGULAR ACCELERATIONS ( R AO I AN S /SEC . S 0 . I 



POOTH 

* 0. 

onniH 

a 

0. 

ANGLES OF 

ATTACK. SIDESLIP (»AnT«NC) 



ALPHAH 

■ .639227F-OI 

R FT AH 

■ 

-.363880E-03 


ZH • -• 139000E+0A 

ZDOTH • 0. XDOTVH • .103000E+03 

WH ■ .6AA573E+01 

DWDOTH * 0. 

PSIH • .975AA3E+02 

RH • -.9289B1E-05 

ROOTH • 0. 


RATE GYRO OUTPUTS 

PI » -.32159F-03 

01 « .2<.612E-02 

PI * .85570F-03 

T IMF • • 1 0000E+02 


LINEAR ACCELFROMETFR OUTPUTS 


UOOTI * .23562E+01 
VOOTI - -.59122E+00 
WOOTI - -.32061E*02 


FEET 

DEGREES 


TINE 


10000E+02 
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LOCAL LEVEL NAVIGATOR 


MISALIGNMENT 

FPR09S 

POSITION 

ERRORS 


VELOCITY ERRORS 



NORTH 

(RAO) 

fast hpwn 

(RAO) (Ran) 

l A T 
(RAD) 

LONG 

(RAD) 

ALT 

(FT) 

LAT RATE 
(RAD/SEC) 

LONG RATE 
(RAD/SEC) 

VERTICAL " 
(FT/SEC) 


-.70694E-03 

.45655E-0* « ft B E-02 

•A2091E-07 .A2910E-07 


I .90960E-08 

.91608E-08 


I ' 


FINAL GUIDANCE COMMANDS 
TIMF » • 10000F+02 SEC 

XDOTGSC (FT/SEC) * .P9500Ft02 PSIHRC (DEG) * .9Q000E+02 ZDOTGSC (FT/SEC) * .988A9E+01 

ZGSC (FT) • -.LB150E+PA 

PHIGSC ( DFG I • 0. PSIGSC (DEG) * .90000E+02 


CONTROLLER OUTPUTS (INCHF«) 

DCC(l) DCC ( 2 ) nrr(3) DCCIA) 

.29038E-02 -.65522E-01 .17A33E+00 -.17019E-01 

TIME - .10000E+02 
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TIME ■ 


20000E+02 


ZH ■ -.131212E+CK FEET 

DEGREES 

ZDDTH » .96<J572Em XDOTVH • .912234E+02 

WH * • 1B7092E+02 

DWDOTH • .1876R7E+00 

PSTH ■ .8973*^+02 

RH ■ .B9?*7lE-03 

RDOTH • -,2865<>3E-03 
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LOCAL LEVEL NAVIGATOR 


MISALIGNMENT 

ERRORS 


POSITION 

ERRORS 


VELOCITY ERRORS 


NORTH 

EAST 

nnwN 

LAT 

LONG 

ALT 

LAT RATE LONG RATE 

VERTICAL 

(RAD) 

(RAD) 

(»*n) 

(RAD) 

(PAD) 

(FT) 

(RAD/SECI (RAO/SEC) 

(FT/SEC) 

-• 89806E-03 

.88 328E-09 

. 5*3^96-02 

. 18 38 0 E- 

•06 • 19259E-06 


I • 19219E-07 • 1 9229E-07 



FINAL GUIDANCE COMMANDS 
TIME * • 20000E+02 SEC 

XDOTGSC (FT/SEC) * .8OOO0F+0P PSIHRC (DEG) * .90000E+02 ZDOTGSC (FT/SEC) » .99927E+01 

ZGSC (FT) * -.12150E+0A 

PHIGSC (DEG) » 0. "FIGSC (DEG) * .90000E+02 


CONTROLLER OUTPUTS (INCHF*) 


DCC(I) DCC ( 2 ) nrr( 3) 

-.O7A28E-01 -*857102-01 *m72Et00 

TIME ■ .20000E+02 


DCC (8 ) 
16618E-01 


UJ 

VO 
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ZH * 121882E + 0<. FEET 

DEGREES 

ZDDTH - .R792R5E*01 XDOTVH « .813297E»02 

WH » • 1R6583E+02 

DWDOTH ■ -• 198721 E-Ol 

PSIH « • 899630E+02 

RH - .119R57E-03 

RDOTH ■ -•381376E-0R 
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